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In this paper we consider quantum metastability in a class of moving potentials introduced by 
Berry and Klein. Potential in this class has its height and width scaled in a specific way so that 
it can be transformed into a stationary one. In deriving the non-decay probability of the system, 
we argue that the appropriate technique to use is the less known method of scattering states. This 
method is illustrated through two examples, namely, a moving delta-potential and a moving barrier 
potential. For expanding potentials, one finds that a small but finite non-decay probability persists 
at large times. Generalization to scaling potentials of arbitrary shape is briefly indicated. 
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I. INTRODUCTION 



An interesting issue in cosmology is the evolution of metastable states in the early universe according to the original 
version and its variants of the inflationary models 0,0- I n these models inflation of the early universe is governed 
by a Higgs field trapped in a metastable state. Inflation ends when the metastable state decays to the true ground 
state of the universe. During inflation the universe expands exponentially. It is thus obvious that the metastable 
state of the Higgs field is trapped in a rapidly varying potential. The problem is therefore a truly time-dependent 
one. However, owing to the inherent difficulties of the problem, more often than not one considers the decay of the 
Higgs field in a quasi-stationary approximation, in which the decay is studied assuming a static potential Q . Surely 
this approximation is hard to justify, but for the present one has to be content with it. Ultimately one hopes to be 
i-H ' able to tackle the non-stationary case. To this end, it is desirable to gain some insights first by studying metastability 
, in time-dependent potential in simple quantum-mechanical models. 

Roughly speaking, time-dependent potentials can be divided into three classes. In the first class we have potentials 
with time-dependent strength. When the strength is small, the Schrodinger equation can be solved by time-dependent 
perturbation theory. Almost all textbook examples belong to this type. When the strength of the potential is not 
' small, other methods of solution must be sought. For example, solutions of time-dependent harmonic oscillator £| 
and time-dependent linear potential |5( can be obtained by the method of invariant. We note here that the interesting 
phenomenon of quantum tunneling induced by an externally driven field has also been examined experimentally and 
theoretically 0,0)0]. The second class of potentials involves time-dependent boundaries. Unlike the first class, this 
class of potentials attracts much less attention, and almost all previous works in this area concerned only the simplest 
of all cases, namely, an infinite potential well with a moving wall 0, The last class is the combination of the 
JL" 1 previous two classes. 

We believe that the barrier potential in an inflationary universe is non-stationary, not only the barrier height but 
also the barrier width should be changing as time elapses. However, it will be extremely difficult to study metasbility 
in such a time-dependent potential in full generality. Thus it would be helpful if the quantum tunneling effect could be 
studied in any class of moving potential, special though it is, as a step to understanding the decay of a non-stationary 
metastable system. 

In this paper we consider quantum metastability in a class of scaling potentials which allows one to apply techniques 
used in the corresponding problem with stationary potentials. This class of potentials was introduced by Berry and 
Klein 11]. Potentials in this class have their heights and widths scaled in a specific way so that one can transform 
the potential into a stationary one. 

The organization of the paper is as follows. In Sec. II. we give a general discussion of the solutions of the Schrodinger 
equation with the scaling form of time-dependent potential introduced in 11]. It is argued that the most suitable 
technique for studying quantum metastability in such kind of potential is the less known method of scattering states. 
Two simple examples of such metastable systems, a moving delta-potential and a moving square barrier, are inves- 
tigated in Sect. Ill and IV, respectively. Generalization to arbitrary barrier is briefly discussed in Sec. V. Sec. IV. 
concludes the paper. 
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II. SCHRODINGER EQUATION WITH A SCALING POTENTIAL 



We shall consider the problem of quantum metastability of a particle of mass m trapped in a moving potential 
V(x,t). We assume that the potential V(x,t) is of the scaling form proposed by Berry and Klein [Tl|, namely, 
V(x,t) = V(x/ L(t))/ L 2 (t), where L(t) is a time-dependent scaling factor. The Schrodinger equation is 



ih 



dV(x,t) 

at 



h 2 d 2 



1 



V 



2m dx 2 L 2 (t) \L(t) 



(1) 



So far solution of Eq.Q is restricted mostly to the special case in which V has the functional form of an infinite 
potential well, i.e. V(x,t) is an infinite well with a moving wall |UEJ- In this case the scaling factor L 2 (t) in front 
of V is immaterial. 

For our purpose, we shall assume V to have the generic shape of a potential well which is impenetrable to the left 
and has a finite barrier to the right, much like that usually employed in the discussion of alpha-decay. We further 
assume L(t) to be a linear function of time: 



L(t) — Lq + vt , v = constant. 



(2) 



Of course, for v < 0, the problem is meaningful only for time duration < t < Lq/\v\. Eq.© cannot be solved 
by separating the time and spatial coordinates. However, for the scaling form of V(x,t) in Eq.© and the linear 
form of L(t), separation of variables can be achieved through a series of transformations introduced in 0,0] (see 
also: jUEl)- O ne first transforms the coordinate frame into a rescaled frame with a rescaled coordinate x defined by 



In this frame the Schrodinger equation becomes 



d 



x(t) 



i) 2 



L(t) ■ 



v _ d - 
i^-rx— + V(x) 
L ox 



2mL 2 dx 2 

Eq. (0J can be further simplified by the following transformation 

= -jL= e^ LvS;2 ^{x,t) 



*(x,t) 



and the introduction of a new time variable r 



ds 



o L 2 {s) L L(t) 
After substituting Eq.© and © into (@J), we obtain the equation 

d Ti 2 d 2 _ - 

ih—$(x, t) = t) + V{x)$(x, t) 



(3) 



(4) 



(5) 



(6) 



(7) 



which resembles the Schrodinger equation with a stationary potential. Eq. (JJJ) can be solved by separation of variables 

$(x,r) = <P{x)e~^ ST , (8) 

where &(x) satisfies the eigenvalue equation 



2m dx 2 



$ fc (x) = E k § k {x) . 



(9) 



Once Eq.© is solved exactly in the rescaled frame, the exact wave function in the original frame is then given by 

1 



*fc(M) 



(10) 
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The set of solutions (|10fl is complete and orthonormal 

($ k {x,t)\*i{x,t)) = = Su , (11) 

so using this set of solutions we can find a solution satisfying any initial condition. Furthermore, if an initial state 
^(x,0) is expressible in the basis {^k} as 

V(x,0)=J2 c k*k(x,0) , c k = {V k (x,0)\V(x,0)) , (12) 
fe 

then at a later time t the state is 

=J2 c ^k(x,t) . (13) 
fe 

We have now succeeded in transforming the original time-dependent Schrodinger equation into a time-independent 
one. The problem of calculating the decay probability of a particle confined in V(x,t) at time t is reduced to the 
corresponding problem with a static potential V(x). Hence techniques used in the time- independent potential for 
calculating decay rate can be borrowed. 

However, there are some subtleties. Naively, one is tempted to employ the most well-known method, namely, the 
complex eigenvalue method, proposed by Gamow in his studies of the alpha-decay [13J. In this approach an "outgoing 
wave boundary condition" is imposed on the solutions of the Schrodinger equation for the particle trapped in the well. 
That means incoming plane wave solutions outside the potential well are discarded right from the beginning. This 
procedure naturally leads to an eigenvalue problem with complex energy eigenvalues. One then relates the imaginary 
parts of the energy to the decay rate. While the c omp lex eigenvalue method is straightforward and physically 
reasonable, it suffers from some conceptual difficulties 14]. For example, how can energy eigenvalues be complex as 
we are dealing with a Hermitian Hamiltonian? Also, the eigenfunctions are not normalizable, a difficulty directly 
related to the eigenvalues being complex. Furthermore, the particle trapped in the well cannot be in an eigenstate of 
the system in the first place, since such states are not completely confined at t = 0. 

Apart from the difficulties mentioned above, the complex eigenvalue method cannot be employed in our problem for 
other reasons. First, the problem we are interested in is an intrinsically time-dependent one, with a non-conservative 
Hamiltonian. Hence energy eigenvalues and eigenstates lose their meanings altogether {E in Eq.JSJ) is not an energy 
eigenvalue). Second, the "outgoing wave boundary condition", essential to Gamow's method, cannot be imposed in 
our case. The reason is as follows. As we discussed before, in order to fix a moving potential we need to transform 
our problem to a corresponding static one in a rescaled frame. But in this frame the meaning of incoming or outgoing 
plane wave is rather obscure. In fact, it can be checked that an outgoing plane wave in the original x-t (rescaled x-t) 
frame contains both "incoming" and "outgoing" components in the rescaled x-t (original x-t) frame. 

Instanton method is another technique commonly used in the calculation of the decay rate of a metastable state 
(lEf . This semiclassical method amounts essentially to finding the imaginary part of the ground state energy of the 
system. Again, it is not suitable for our case for the same first reason given above for the failure of the complex 
eigenvalue method. 

Now that the two most common methods fail to suit our purpose, we have to look for alternatives. Fortunately, a 
different method exists, namely, the scattering state method (or virtual level method, as Fermi called it) [16J. This 
method is much less well known and seldom used in the literature 0. However, it is conceptually the most satisfying 
one of all the methods. In this method, one first constructs the initial confining state, which is not viewed as an 
eigenstate, as a linear superposition of scattering states with real energies, and follows its evolution in time. In the 
course of this evolution, no energy will become complex. Unlike the Gamow states, the scattering states contain both 
incoming and outgoing components in the region into which the particle escapes. It is this feature of the method 
that makes it most suitable for our present problem. The method is easily adapted to Eq. Q) by taking the scattering 
states as the states © with real values of E. 

In the following two sections, we apply the scattering state method to two simple examples of the class of scaling 
potentials discussed in this section. As the scattering state method is not so well known in the literature, we think it 
appropriate to give some details in order to make our paper self-contained. We shall follows the procedures given in 
|14| which are slightly adapted to our needs. 



III. MOVING DELTA-FUNCTION POTENTIAL 



Our first example is a uniformly moving delta-function potential 



V(x,t) 



00 

,(t) 



r%5(x-a(t)) 



a) 



x < ; 
x > , 



(14) 



4 



where a(t) = aL(t) > gives the location of the delta-potential. This is of the class of potentials defined in the last 
section and corresponds in the rescaled frame to V(x) = oo for x < and V(x) = VqS(x — a). We must have = 
in the region x < 0. In the region x > 0, Eq.Q is 



h 2 d 2 $(x) 
2m dx 1 



V S(x - a)<Z>(x) = E$(x) 



Its general solutions are 



J sin(fca:) 



Ccos(kx + 9) 



< x < a ; 
a < x , 



(15) 



(16) 



where k = y2rnE/h, C is a real constant, and 9 a phase angle. Note that we have chosen the wave function to be 
real, and included incoming wave component in the region x, > a. This ensures that k, and hence E, is always real. 
The wave function and its first derivative satisfy the following boundary conditions at x — a: 



<5>(x 



") = §(x = a ) 



and 



dx 



dx 



2m - 

—2-V $(x = a) 



(17) 



(18) 



^From these relations the coefficient C can be determined as a function of k, 



C z (k) — sin^fca) + ( cos(ka) H sin(fca) 

h ka 



(19) 



Physically, the value of C 2 (k) can be interpreted as the ratio of the probability of finding particles in the region x > a 
to the probability within the confined region < X < a for a particular k. The general shape of C 2 (k) is shown in 
Fig. 1, from which we can assert that the particle can be trapped within the confined region only when C 2 assumes 
one of its minima, which occur only in the neighborhood of some specific values of k. In these regions the values of 
C 2 (k) are extremely small. From Eq. l|19H it is obvious that these minima will center around k n = nir/d (n = 1, 2, . . .) 
(i.e. sin(fc„a) = 0) as long as Vb is large enough so that 2m Vo/h 2 » mr. 

Below we shall restrict our discussions to the case of large Vb ■ In this case approximate analytic expressions can be 
obtained and compared with the corresponding results in the time-independent case [Tsj . According to the scattering 
state method, one constructs confining states in the potential well by taking suitable superposition of the scattering 
states with k in the neighborhood of k n . To this end, let us first expand C 2 about E n — h 2 k 2 l /2m (we revert to the 
variable E below) 



C 2 (E) 



2mV 



(E-E n + sy 



i+feBY 

V h 2 k n J 



G 2 (A + S) 2 + F 2 , 



where A = E — E n and the constants 



2Vo 
a 



G 



F 



2mV 
Tt k n 



2mV 



(20) 
(21) 

(22) 
(23) 
(24) 



The scattering states with E in the neighborhood of E n can then be written as 



■0a(z) 



y/G 2 (A+8f+F 2 



-| cos(fca: + 9) 



sin(fca;) , < x < a ; 



(25) 



a < x . 
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We have quantized our system in the interval [0, R], where R^> a which at the end of the calculation will be set to 
infinity. With these scattering states we now construct an initial state which is completely confined within the well 
by taking linear combination of the scattering states with different A but the same value of E n , 



$(x, t = 0) = ^2 c a^a(x) = 



t (x) , x < a 







x > a . 



The coefficient ca can be calculated from orthogonality of the states -0A (x) , 



Ca 



dxipA{x)$(x,0) 



Choosing 



we get 



4>n{x) ~ \/ -= SHI 



f 



R^G 2 (A + 5) 2 + F 2 Jo 

. /rvKX\ „ „ „ 

im ^ _ J ; n = 1,2,3,... 

dx sin (/ex) sin ( 



dx sm(kx)(j) n (x) 



'RE y/G 2 {A + S) 2 + F 2 Jo 

a I 
Ry/G 2 {A + 5) 2 + F 2 ' 



The initial state is then given by 



4> 



R a V G2 ( A + S ) 2 + F 2 



iPa{x) . 



^From Ea. HlU|) . H12|) and (|I3[) . the solution at a later time r is 

a v - f 



ipA(x)e 



-|(E„+A)t 



(26) 

(27) 

(28) 

(29) 
(30) 

(31) 

(32) 



As the system is quantized in the interval [0, R], we have kR = n'n/2 where n' is a very large integer (n' 3> n). After 
replacing the sum by an integral 



E 

A 



Ea.i|52*)l becomes 



&(x, r) 



R 12m 
nh V E„ V R 



dA 



R F2m 
dA —\-=r 
Trn V E n 



1 



ipA{x)e >> 



-±{E n +A)n 



(33) 



(34) 



V / G 2 (A + S) 2 +F 2 

Substituting V'a(S) from Ea. i|25|) . we obtain the approximate wave function of the state confined in the well (0 < x < a) 

as 



$(je,t) ~ \\hgr sm(k n x)e * E " T 



dA 



G 2 (A + 8) 2 + F 2 



2 I ma sin(fc r) a:) < s „_ iif 



(35) 



For metastable system an important quantity is the non-decay probability P(t) that the particle is still in the well 
at time t if it is initially confined in the well at t — (P(t = 0) = f ). In our case P(t) is defined as 



f a{t) Mx,t)\ 2 dx 
f a{0) \*(x,0)\ 2 dx 
J Q a mx,r)\ 2 dx 
J a mx,0)\ 2 dx ' 



P(t) 



(36) 
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/From Eq.J25J) we find 



P(t) ~ exp(- 7 „(t)) 



where 



7n(*) = 2 



LoL(t) 



2mVb 

7l fcrj. 



L L(i) ■ 

We have used Eq.@, (j2"5)l and (2^ to obtain the result We note that for expanding potential (i> > 0) 



7»(*) 



1 

L w 



(37) 

(38) 
(39) 

(40) 



as t — > oo. Unlike the stationary case (u = 0), there is a small but finite probability that the particle does not tunnel 
out of the well. This result is reasonable, since as the barrier moves away from x = it leaves more room for the 
particle to stay within the well. 

For Vq much larger than the characteristic value of k n of the escaping particles, j(t) becomes 



7n(t) 



(h 2 k n f 

2m 3 aU 2 



t 



L L(t) 



h 6 {nir) 3 
2m 3 a 4 U 2 



t 



L L(t) ' 



(41) 



where k n = nir/a has been substituted. It is proper to compare Eq. I|41|l with the corresponding result in the 
stationary case (v = 0). In the limit v — > 0, we have L(t) — > Lo, a — > a/Lo, E n — > E n Lg, and £;„ — > k n Lo, where E n 
is the corresponding energy in the static frame, and k n = \/2mE n /h. In this limit j(v = 0) is directly proportional 
to the time t. We can therefore define a decay rate by T n = j n (v = 0)/t, which in this case is 



61.3 



h 6 {nir) s 



2m 3 a(^) 2 2m 3 a 4 (^) 2 



(42) 



Ea. (|42|l is consistent with the result obtained by the complex eigenvalue method in 01 f° r a static delta-potential 
located at x = a with strength Vq/Lq. 



IV. MOVING SQUARE BARRIER POTENTIAL 

Next, we consider a moving barrier potential 

{00 , x < ; 

j^Vq , a{t) <x< b(t) ; (43) 
, x > b(t) , 

with a(t) = aL(t) and b(t) = bL(t) (a and b are two positive constants). When the problem is transformed to the 
rescaled frame, it is equivalent to solving Eq. @ with a stationary potential 




V(x) = < U , a < x < b ; (44) 



The general solutions are 



sin(fcx) , < x < a , k - 

= { Ae k'* + Be -W* a <x<b, k> = < / ^(Vo-E) . i U, , 



2mE 



C cos{kx + 8) , a < x 



h 
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Here A, B and C are real constants, and 9 a phase angle. As before, we have set the solutions real in the whole region 
to ensure that k is always real. We also require the solutions and their derivatives be continuous at the boundaries 
x — a and b. These boundary conditions determine the values of the coefficients A, B and C as functions of k: 



A(k) 



1 

r 



-k a 



sin(fca) + -=r- cos(fca) 



B(k) = \e k '-* 



sin(fca) — — cos(fca) 



h' 2 



2k'b \2 



k 2 







) AB + 









,—2k b o2 



(46) 
(47) 
(48) 



The general shapes of A(k) and C 2 (k) are shown in the Fig. 2. We see that metastable states of the system will occur 
only in a finite number of neighborhood of k n (n = 1, 2, . . .) such that A{k n ) = 0. The roots k n satisfy 



- k - 
sin(fc„a) + ■=— cos(fc„a) = . 

If' 



(49) 



Ea. (|48t implies that C 2 (k) is minimal at k n . For a given Vq, the number of roots k n is restricted by the condition 
that k' n in Ea. l{4*5l) . 



'2mV 



k3. 



(50) 



must be real. Hence the possible values of k n can only lie in the interval (0, y 2 ™¥" ) ■ For instance, there are only two 
roots for the parameters assumed in Fig. 2. 

Let us now expand the coefficients A(k) and B(k) about E n (= h 2 k 2 /2m) 



^ - [i 

B{E) » B(E n ) 



(E - E n ) 



E=E n 



(51) 
(52) 



Inserting Ea. f5T|l and l(5^|l into Eq. |@H|l ; an d after some tedious calculations we find that C 2 (E) behaves in the 
neighborhood of E n as 



C 2 {E) = G 2 {A + 5) 2 +F 2 , 
where A = E — E n as in the previous example, and the constants in the present case are 



(53) 



G 



F 



1 / ma 



4 \h 2 kn 



1J2 



1 



k 2 



cos(/c n a) — sin(fe n a) 



a 2k' n (b-a) . 



sin(fc„a) — =^ cos(fc„a) 



-2k' n (b-a) 



and 



h 2 k r 



k 2 - V 2 
h 2 4- Tc' 2 



k' n sin(fc„q) - k n cos(fc„q) ^ p _ 2 fc^(b-a) 
k' n cos(fc„a) - fc„ sin(fc„a) / 



(54) 
(55) 

(56) 



Once we obtain the relation (|53|l . we can construct the scattering states relevant to this metastable system following 
exactly the same procedures as those in the previous section. The non-decay probability P(t) of finding the particle 
within the confined region (0 < x < a) at time t is again of the from P{t) ~ exp(— 7 n (t)), where 7 n (t) is now given by 



7n(*) = 2 



F 
G 

8h 2 kl 



t 



L L(t) 



L.2 _i_ L/2 



-2fe^(fc-a)_ 



L>oL(t) ' 



(57) 
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where we have used Eq . (|54|) . i|55|) and (|49|l . When the barrier height is much larger than the characteristic "energy" 
of the escaping particles (Vq E n ), which is equivalent to the relation k' n 3> k n , j n (t) becomes 

ln(t) — ^§ e -^(5- a ) * . (58) 
ma /c„ L^Lyt) 

As in the previous example, for positive v > (the expanding case) one finds a small but finite probability that the 
particle does not tunnel out of the well at large time . In this case not only does the barrier leave more room for 
the particle to stay within the well as it moves away from the x = 0, but its width also become thicker, thus making 
tunneling difficult. 

In order to transform the result (|58|l to the stationary one, we again make use of the same substitutions as given at 
the end of the last section, with the addition of b — > b/L and k' n — » k' n L , where k' n = y/2m((Vo/Lo) — E n )/h. Once 
again 7„(w = 0) is directly proportional to the time t, in which case a decay rate can be defined: T n = 7 n (« = 0)/t. 
For the present example we have 

Tn (v = 0) = — ^L e -2K(b-*) (59) 
ma k'^ 

which is the same as the result obtained by the complex eigenvalue method for a square barrier with width (b — a) 
and height PS- 



V. GENERAL SCALING POTENTIALS 



We have calculated the non-decay probabilities of two non-stationary metastable systems explicitly. The potential 
barriers in the rescaled frame considered in these systems assumed the form of a delta-function and a square barrier. 
These calculations can be immediately generalized to barriers with more general shapes. Without giving further 
examples, what we would like to do here is to discuss briefly a close connection between the the non-decay probability 
Pit) of a particle in a metastable scaling potential V(x,t) — V \x / Lit)) j 'L 2 it) and the decay rate V of the same 
particle if it were instead confined in a static potential well V(x) = Vix). 

^From the discussions and examples in the previous sections, we know that the calculations of P(t) is reduced to the 
corresponding computations in a static potential Vix) in the rescaled frame. Now the later task would be exactly the 
same as that carried out in the potential Vix) = Vix) in ordinary coordinates. The only difference, as seen from the 
previous two examples, is that all ordinary parameters, such as E, k, k', t, a, etc, are replaced by the corresponding 
rescaled ones, i.e. E, k, k', r, a, etc. Application of the scattering state method to the general alpha-decay type 
of potential Vix) in normal coordinates has been given in [T^j . and can be carried over directly. Following the 
important step is to determine the discrete values E n (or equivalently k n ) that minimize the amplitude C of the wave 
function in the region outside the well. Consider a confining state constructed with E centered around a specific E n . 
Minimization of C then gives the two functions FiE n ) and GiE n ) (other parameters in F and G are not indicated). 
The non-decay probability in Vix) is then given by exp(— T n t), where the decay rate r„ is 



^n{E n ) — 2 



F(E n ) 



GiE n ) 



(60) 



Suppose all these computations have been done in ordinary coordinates. Then one can immediately write down the 
expression of the non-decay probability P(t) ~ exp(— 7„(i)) for the scaling potential Vix,t) as 



7n(*) 



F(E n 



GiE n ) 

r„(£n) 



L L(t) 



L L(t) 



(61) 



Here the functional form of the decay rate T n is taken over directly, but with all the parameters replaced by the corre- 
sponding rescaled ones. Eq.JHH gives the connection between the non-decay probability in Vix, t) = Vix/ Lit))/ 1? it) 
and the decay rate in Vix) = Vix). Finally, we note here that, in the non-moving limit v = 0, Vix, t) becomes 
Vix) = V ix / Lq) / L/q. Setting v = in Eq. ljBTTl then gives the decay rate in this potential: r n (P n )/Lg, as we had seen 
in the previous cases. 
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VI. CONCLUSION 



In this paper we consider quantum metastability in a class of moving potentials introduced by Berry and Klein. 
Potential in this class has its height and width scaled in a specific way so that it can be transformed into a stationary 
one. In deriving the non-decay probability of the system, we employed a method which is less well known but 
conceptually more satisfactory, namely, the method of scattering states. Non-decay probabilities in a moving delta- 
potential and a moving square barrier potential were derived. We also give a connection between the non-decay 
probability in a general scaling potential and the decay rate in a related static potential. 
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Figure Captions 

Figure 1. The shape of In C 2 {kd) in Ea. l|19[) as a function of ka for 2mVo/h 2 = 10 (dotted line) and 200 (solid line). 
The minima of In C 2 (ka) will center around k n a = wk (n = 1, 2, . . .) for large values of 2mVo/H 2 . 

Figure 2. The shapes of In C 2 (ka) (solid line) in Eq. ljlgj) and 5e^' a A(ka)/2 (dotted line) in Ea. (|^S|l as function of 
ka for b = 2d and 2mVoa 2 /h 2 — 40 . It shows that the minima of In C 2 (ka) only occur in a finite number of 
neighborhood of k n (n = 1, 2, . . .) such that A(k n d) = 0. 



